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A CLT CONCERNING CRITICAL POINTS OF RANDOM FUNCTIONS 

ON A EUCLIDEAN SPACE 

LIVIU I. NICOLAESCU 


Abstract. We prove a central limit theorem concerning the number of critical points in 
large cubes of an isotropic Gaussian random function on a Euclidean space. 
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1. Introduction 

Throughout this paper X(t) denotes a centered, isotropic Gaussian random function on 
m > 2. 

Assume that X is a.s. C^. For any Borel subset S C we denote by Z{S) the number 
of critical points of X in S. For a positive number L we set 

ZL:=Z{[-L,Lr), 

and we form the new random variable 

( 1 . 1 ) 
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In this paper we will prove that, under certain assumptions on X{t), the sequence of random 
variables (Cn), converges in distribution as —)• oo to a Gaussian with mean zero and finite, 
positive variance. 

The proof, inspired from the recent work of Estrade-Leon [13], uses the Wiener chaos 
decomposition oi Z^. 


Notations. 


• N Z>o, No Z>o. 

• For any positive integer n we denote by the identity map M” —>• M”. 

• For t = {ti,..., tm) £ we set 


t 


m 



max \tk\- 
l<k<m 


• For any u > 0 we denote by 7 ^ the Gaussian measure on M with mean 0 and variance 

V. 

• If X is a scalar random variable, then we will use the notation X £ N(0, v) to 
indicate that X is a normal random variable with mean zero and variance v. 

• If A is a subset of a given set 5, then we denote by I a the indicator function of A 


Ia:S^ {0,1}, Ia{s) 


1, s £ A, 

0, seS\A. 


• If X is a random vector, then we denote by E[X] and respectively var{X) the mean 
and respectively the variance of X. 


2. Statement of the main result 
Denote by K{t, s) the covariance kernel of X{t), 

K{s,t) := E[X{t)X{s)], t,s£M™. 

The isotropy of X implies that there exists a radially symmetric function C : M”* —)• M such 
that K{t,s) = C{t — s), yt,s £ M. We denote by fi{dX) the spectral measure of X so that 
C{t) is the Fourier transform of // 

C'(t) = (27r)-f [ (2.1) 

JR'" 

2.1. The setup. For the claimed central limit result to hold, we need to make certain as¬ 
sumptions on the random function X{t). These assumptions closely mirror the assumptions 
in [13]. 

Assumption Al. The random function X{t) is almost surely C^. 

To formulate our next assumption we set 

' ip { t ) := msix {\ dfC { t )\-, ]al<4}, t £ M”*, (2.2) 

where for any multi-index a = (oi,..., am) £ N™ we set 

|«| :=ai + ... + a^, df = . 


lim = 0 and -0 £ L^(M™'). 

|t|^oo 


Assumption A2. 
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Our next assumption involves the spectral measure /x(dA) and it states in precise terms that 
this measure has a continuous density that decays rapidly at oo. 

Assumption A3. There exists a nontrivial even, continuous function re : M —)■ [0, oo) such 
that 

//(dA) = u;(|A|)dA. 

Moreover w has a fast decay at oo, i.e., 

|A|^u;(|A|) e 

Remark 2.1. (a) Let us observe that A1-A3 imply that 

-if G \/q > 0. 

(b) The assumptions A1-A3 are automatically satisfied if the density w is a Schwartz function 
on M. 

(c) The paper [13] includes one extra assumption on X, namely that the Gaussian vector 

J2(A(0)) := (A(0),VA(0),V2a(0)). 

is nondegenerate. We do not need this nondegeneracy in this paper, but we want to mention 
that it is implied by Assumption A3; see Proposition A. 6 . □ 

Fix real numbers u > 0 and u > 0. Denote by Sm the space of real symmetric m x 
m matrices, and by the space Sm equipped with the centered Gaussian measure T^^^ 
uniquely determined by the covariance equalities 

E[aijak^] — aSijSki T vi^Si^dj^ T VI ^ i-ij^ .k,£ ^ m. 

In particular we have 

E[a‘fj] = u + 2v, E[aiiajj] = u, E[a^j] = v, VI < i 7 ^ j < m, (2.3) 

while all other covariances are trivial. The ensemble § 1 ^ is a rescaled version of the Gaussian 
Orthogonal Ensemble (GOE) and we will refer to it as GOE])^. As explained in [22, 23], 
the Gaussian measures are invariant with respect to the natural action of 0{m) on S^. 
Moreover 

dro,^(A) = e-^^^^^\dA\. (2.4) 

The ensemble can be given an alternate description. More precisely a random A G 
can be described as a sum 

A = B + Aim, B G GOEm, X G N{0,u), B and X independent. 

We write this 

Sr = GOE];+Ar(0,u)lm, (2.5) 

where + indicates a sum of independent variables. We set Sm := Recall from (2.1) that 

E;[ A(t)A(0) ] = C{t) = K{t, 0) = (27r)-t [ e-^^*’^'^w{\X\)dX. 

Eollowing [25] we define 

Sm ■= ;- TK I w(\x\)dx, dm ■= ;- 777 / x'^w(\x\)dx, 

(27r)-/2y^^ VI II , m (27r)™/2yjj^ 1 

hm ■= ;- 777 / x‘ixnw(\x\)dx. 

(271)^/^ 1 2 VI II 
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Clearly Sm, dm, hm > 0. If we set 

poo 

Ik{w) := / w{r)r^dr, (2.6) 

Jo 

then we have (see [25]) 

m m 

" (2.7) 

Then we deduce that 


t;[x(o)-4X(o)] =T;[4X(o)-a2.,^x(o)] =o, Vz,j,/c ( 2 . 8 a) 

E[X{0f] = Sm, E[duX{Q) • 5t^X(0)] = Vi,i, (2.8b) 

E[X{Q) • 52i^.X(0)] = Vz,j, (2.8c) 

T;[ t.X(0) • dl,^X{Q) ]=hm{ Sijdke + SikSji + SieSjk), Vi, j, k, L (2.8d) 

The equality (2.8b) shows that VX(0) is a M^-valued cenetered Gaussian random vector 
with covariance matrix dm^m, while (2.8d) shows that V^X(O) € □ 

2.2. The main result. We can now state the main result of this paper. 


Theorem 2.2. Suppose that X{t) is a centered, stationary, isotropic random function on 
M™’, m > 2 satisfying assumptions Al, A2, A3. Denote by Z^- the number of critical points 
of X{t) in the cube Cj^ := [—A, A]™. Then the following hold. 

(i) 

E[ZM]=Cm{w){2Nr, '^N, (2.9) 

where 

(ii) There exists a constant Voo = Voo{m,w) > 0 such that 

var{ZN) ~ IVogA™ os A —)■ oo. (2-11) 

Moreover, the sequence of random variables 

Cn = N-^/^(^Zn-E[Zn]) 

converges in law to a normal random variable Coo with mean zero and positive vari¬ 
ance Voo- 

Remark 2.3. (a) The isotropy condition on X{t) may be a bit restrictive, but we believe 
that the techniques in [13] and this paper extend to the more general case of stationary 
random functions. (The only place where isotropy plays an crucial role is in the proof of 
Proposition 3.1.) However, for the geometric applications we have in mind, the isotropy is a 
natural assumption. Let us elaborate on this point. 

Suppose that (M, g) is a compact m-dimensional Riemannian manifold, such that volg(M) = 
1. Denote by p the injectivity radius of g. For e > 0 we denote by pe the rescaled metric 
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:= e ‘^g. Intuitively, as e ^ 0, the metric becomes flatter and flatter. Denote by 
the Laplacian of g and by of g,.- Let 

^0 ^ -^1 ^ -^2 ^ • 

be the eigenvalues of Ag, multiplicities included. Fix an orthonormal basis of L?‘{M,dVg) 
consisting of eigenfunctions 'Lfc of Ag, 

Ag^k = 

Then the eigenvalues of Ag are Afc(e) = £^\k with corresponding eigenfunctions 'Ll := 

We now define the random function Y^{p) on M, 

,_ 1 

Yeip) = Xk{e))^ Zk^lip), 

k=0 

where {Zk)k>o is a sequence of independent standard normal random variables. 

Fix a point Pq £ M and denote by exp^ the exponential map Tp^M —)• M defined by the 
metric g^. (This is a diffeomorphism onto when restricted to the ball of radius e~^p of the 
tangent space Tp^M equipped with the metric g^.) Denote by IRg the rescaling map 

Tp^M Tp^M, u i-A eu. 

This map is an isometry {Tp^M,g) —>• [Tp^M, g^). We denote by the random function 
on the Euclidean space {Tp^M,g) obtained by pulling back via the map exp^oIRg. The 
random function X^{t) is Gaussian and its covariance kernel converges in the C°° topology 
as e —)• 0 to the covariance kernel of random function X we are investigating in this paper. 
Denote by N{Xe, Bj.) the number critical points of in a g-hal\ of radius r < p on M, and 
by N{X, Br) the number of critical points of A in a ball of radius R. In [25] we have shown 
that 

E[N{Ye,Br)] ^ E[N{X,B^/,] = const. £-^ as e ^ 0. 

Additionally, in [24] we looked at the special case when M is a flat m-dimensional torus and 
we showed that the variances random variables 

g-m/2 (^NiY„Br)-E[ NiY„ Br)]), ( N{X, B,/,) - E [ N{X, B,,,) ] ) 

have the identical finite limits as e > 0. In [24] we were not able to prove that this common 
limit is nonzero, but Theorem 2.2 shows this to be the case. 

These facts suggest that the random variable N(Y^, Bj.) may satisfy a central limit theorem 
of the type proved in [6, 15]. We will pursue this line of investigation elsewhere. □ 

2.3. Outline of the proof. The strategy of proof owes a great deal to [13]. The Gaussian 
random variables X{t), t £ M™, are defined on the same probability space (D,T, P). We 
denote by .SX the L^-closure of the span of these variables. Thus, tZ' is Gaussian Hilbert 
space in the sense of [16, 19], and we have a white-noise Hilbert space isomorphism 

W : 3C. 

We set 3C := P), where T C T is the cr-algebra generated by the variables X{t), 

t £ M™. __ 

Using the Kac-Rice formula we show in Subsections 3.1 and 3.2 that Z^^ £ ^ and we 
describe its Wiener chaos decomposition. In particular, this leads to a proof of (2.9). To 
prove that the random variables 


Cn = N-^/\Zn - E[Zn]) 
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converge in law to a normal random variable we rely on the very general Breuer-Major type 
central limit theorem [26, Thm. 6.3.1], [27]. Here are the details. 

Denote by the g-th chaos component of (jy G ^. According to [26, Thm.6.3.1], to 
prove that (ly converges in law as —)■ oo to a normal random variable Coo it suffices to 
prove the following. 

(i) For every g G N there exists Vq,oo > 0 such that 

lim var{CN,q) = Vq^oo- 
N^oo 

(ii) 

lim sup y var{(^]y„) = 0. 

Q^oo jy 

q>Q 

(hi) For each g G N, the random variables CAr,g converge in law as iV —)• oo to a normal 
random variable, necessarily of variance Vq^oo- 

We prove (i) and (ii) in the first half of Subsection 3.3; see (3.13) and respectively Lemma 

3.4. In the second half of this subsection we prove that V 2 ,oo > 0. 

To prove (iii) we rely on the fourth-moment theorem [26, Thm. 5.2.7], [29]. The details 
are identical to the ones employed in the proof of [13, Prop. 2.4]. 

The variance of the limiting normal random variable N^o is 

Poo — ^ ^ pQiOO ^ 
q>2 

Appendix A contains estimates of the lower order terms in the Hermite polynomial decom¬ 
position of j det A] where A G m. S> 1. These estimates can be used to produce explicit 
positive lower bounds for Voo for large m. 

2.4. Related work. Central limit theorems concerning counts of zeros of random functions 
go back a while, e.g. T. Malevich [20] (1969) and J. Cuzik [10] (1976). 

The usage of Wiener chaos decompositions and of Breuer-Major type results in proving 
such central limit theorems is more recent, late 80s early 90s and we want to mention here 
the pioneering contributions of Chambers and Slud [9], Slud [30, 31], Kratz and Leon [17] 
and Sodin-Tsirelson [32]. 

This topic was further elaborated by Kratz and Leon in [18] where they also proved a 
central limit theorem concerning the length of the zero set of a random function of two 
variables. We refer to [7] for a particularly nice presentation of these developments. 

Azais and Leon [6] used the technique of Wiener chaos decomposition to give a shorter and 
more conceptual proof to a central limit theorem due to Granville and Wigman [15] concerning 
the number of zeros of random trigonometric polynomials of large degree. This technique 
was then succesfully used by Azais, Dalmao and Leon in [5] to prove a CLT concerning the 
number of zeros of Gaussian even trigonometric polynomials and by Dalmao in [11] to prove 
a GLT concerning the number of zeros of usual polynomials in the Shub-Smale probabilistic 
ensemble. Recently, Adler and Naitzat [1] used Hermite decompositions to prove a GLT 
concerning Euler integrals of random functions. 

Acknowledgments. I want to thank Yan Fyodorov for sharing with me the tricks in Lemma 

A.5. 
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3. Proof of the main result 


The random variables X{t), t £ M"* are defined on a common probability space (hi, 0, P). 
We denote by Ox the u-subalgebra of 0 generated by the variables X{t), t € M™. For 
simplicity we set := L‘^{Q,0x,P)- 

As detailed in e.g. [16, 19, 28], the space admits an orthogonal Wiener chaos 

decomposition 

OO 

L\n) = ^L^n)„ 

q=0 

where L‘^[Q)q denotes the ( 7 -th chaos component. We let Tq : L^(n) —)• L‘^[Q)q denote the 
orthogonal projection on the g-th chaos component. 

Let T denote a compact parallelipiped 

T := [ai,bi] X ■ ■ ■ X [am,bm], ai<bi, Vi = l,...,m. 

From [2, Thm.11.3.1], we deduce that X is a.s. a Morse function on T. In particular, almost 
surely there are no critical points on the boundary of T. 


3.1. Chaos decompositions of functionals of random symmetric matrices. The dual 
space (Sm)* = Hom(§]jj,M) is a finite dimensional Gaussian linear space in the sense of [16] 
spanned by the entries of a random matrix A £ Its Fock space is the space 

L?‘{$rn-,^v,v) and admits an orthogonal chaos decomposition. 


q>0 


We recall that 



<? 




fc=o 


is the subspace of spanned by polynomials of degree < n in the entries of A £ 

and L?‘{$'^)q is the orthogonal complement of '^q-i^v in The summand L‘^{§'!^)q is called 
the q-th chaos component of ). 

The chaos decomposition construction is equivariant with respect to the action of 0{m) 
on In particular, the chaos components L‘^{§^)k are 0(m)-invariant subspaces. If we 
denote by the subspace of consisting of 0 (m)-invariant functions, then we 

obtain an orthogonal decomposition 


= (3.1) 

k>0 

where consists of the subspace of where 0{m) acts trivially. In particular, 

we deduce that L'^{§^y™ consists of polynomials in the variables tr A,tr A^,... ,tr A”*. 

We define the 0(m)-invariant functions 

p, g, / : SmK, p(A) = (tr A)^, g(A):=trA^, /(A) = |detA|. 

A basis of T 2 ™ is given by the polynomials 1, tr A, p(A), g(A). Clearly, since tr A is an odd 
function of A, it is orthogonal to the even polynomials l,p(A),g(A). We have 

n 

Eg-^ [p(A) ] = E[a‘^j] + 2 E[aiiajj] = ‘imv + m{m — l)v = m{m + 2)v. 

i=l i<j 


(3.2) 





LIVIU I. NICOLAESCU 


We have 

m 

[q{A)] = E[al] + = 3m + m(m - 1) = m(m + 2)v. (3.3) 

2=1 i<j 

We deduce that the polynomials 

p{A) = p{A) - Esv \p{A)] = p{A) - m{m + 2)u, 

13.4) 

q{A) = q{A) - E§v^[qA)] = q{A) - m{m + 2)v 
form a (non-orthonormal) basis of Lp'{$'^'^ 2 '’ ■ 

3.2. Hermite decomposition of Z{T). For e > 0 define 

4 : ^ M, de = 

Observe that the family ((5^) approximates the Dirac distribution 6q on as e \ 0. We 
recall [13, Prop. 1.2] which applies with no change to the setup in this paper. 

Proposition 3.1 (Estrade-Leon). The random variable Z{T) belongs to Lp{Q). Moreover 

Z{T) = lim [ |detV2w(t) |4(VX(t))dt 
e\0 Jt 

almost surely and in L^(D). □ 


The above nontrivial result implies that the random variable Z(T) has finite variance and 
admits a chaos decomposition as elaborated for example in [16, 19, 28]. 

Recall that an orthogonal basis of L^(M, 7 ^(dx)) is given by the Hermite polynomials, 
[16, Ex. 3.18], [21, V.1.3], 


Hn{x) := 


In particular 


Eor every multi-index ol = 
and any 


/ -.-.n — { 

(- 1 ) e 2 -— 

^ ^ dx^\ 


e - =ni 


Ifl 

E 

r=0 


(-ir 


n—2r 


2 ’'r!(n — 2r)!' 


(3.5) 


Jo, n = 1 mod 2, 


(3.6) 


(oi, 02 , ■ ■ ■) G Nq such that all but finitely many o^-s are nonzero, 
x = (xi, X2,...) G 


we set 

\a\:=Y(^k, a!:=]Jafc!, Ha{x) := Y\_Ha^{xk). 

k k k 

To simplify the notation we set 

U{t) := ^VX(t), A{t) := V^X{t). 

V 

Thus U (t) is a M™-valued Gaussian random vector with covariance matrix Im while A{t) is 
a Gaussian random symmetric matrix in the ensemble 

Recall that f{A) = \ detHj. We have / € and we denote by /n(H) the com¬ 

ponent of f{A) in the n-th chaos summand of the chaos decomposition (3.1). Since / is an 
even function we deduce that fn{A) = 0 for odd n. Note also that 

/o(H) = E;g^^[ldetHl] /O. 
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Following [13, Eq.(5)] we define for every a G N™ the quantity 

d(a) := ■^(27rdm)~^ffa(0). (3.7) 

al 

Arguing exactly as in the proof of [13, Prop. 1.3] we deduce the following result. 
Proposition 3.2. The following expansion holds in Lf‘{Tt) 

OO 

Z{T) = Y,Z,{T), 

g=0 

where 

Zg{T) = PgZ{T)= d{a) f H^{U{t))fn{A{t))dt. □ 

aeNg",neNo, 

\a\+n=q 

Observe that the expected number of critical points of A on T is given by 

E[Z{T)] =E[Z{T)]=d{0) [ E[Ho{U{t))fo{A{t))]dt 

Jt 

(use the stationarity of X{t)) 

= {2-Kdm)~^ j E^^[f{A{t))]dt = {27rdm)~^ E^,^[\detA\]vol{T) 

This proves (2.9) and (2.10). 

3.3. Asymptotic variance of Z{T). Denote by Cn the cube [—A, and by B the cube 

[0,1]™'. For any Borel measurable subset S C such that vol(S') 7 ^ 0 we set 

a/vo1(5) 

Thus, (n = C(CAf)- Since Cn € T^(D) we deduce var{fN) < 00 . 

Proposition 3.3. There exists Poo £ (0, 00 ) such that 

lim var{C,N) = Poo- 

N—^OO 

Proof. To prove that the above limit exists we follow closely the proof of [13, Prop. 2.1]. We 
have 

Cn = aCN) = {2N)-^Y.Zg{CN), 

g>i 

Pjv:=^ar(C(Cjv)) = V (2A)-”^F;[^.(CAr)'] . 

— ■^q,N 

To estimate Vg^jsf we write 

Zq{T) = J^Pq{t)dt, 

Pq{t)= Y. d{a)H^{U{t))U{A{t)). 

nENo 
\a\-\-n=q 


where 


(3.8) 
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Then 


Vg,N = {2N)-^ [ E[pg{s)Pg{t)]dsdt 

J Cj^x C]\f 

(use the stationarity of X{t)) 

= {2N)-^ [ E[pq{0)pq{t-S)]dsdt 

J CnX C]\f 

= £:[p,(0V,(«)] p (l - du. 

The last equality is obtained by integrating along the fibers of the map 


Cj\i X C 3 (s, I— t — s G C 2 N• 

To estimate the last integral, we fix an orthonormal basis of the Gaussian Hilbert 

space Hom(S((’",M). We denote by B the vector {hij)i<i<j<m, by A the vector (ajj)i<j<j<m 

both viewed as column vectors of size 

. m(m + l) 

v{m) := dimhm =- - -. 

There exists a nondegenerate deterministic matrix A of size ^[m) x relating A and B, 

A = KB. We now have an orthogonal decomposition 

fn{A) = Y. 

\P\=n 

Let us set 

3,„ := (N^) X (NS'™’). 

We deduce 

pS)= ^ d{a)c{(3)HaiU{t))c{l3)Hp{B(t)). 

We can further simplify this formula if we introduce the vector 

Y{t):={U{t),B{t)), B{t) = K-^A{t). 

For -y = (a, fd) G we set 

a( 7 ) := d(a)c(/3), H^{Y{t)) := H^{U{t))Hp{B{t)). 

Then 

Pq{'t)= Y “( 7)^7 (^(*))> ( 3 - 9 ) 

l&m, \l\=q 

E[p,{^)p,{u)]= Y <i)<i')E[H,{Y{Q))H,,{Y{u))]. 

J,Y&m 

hl=h'l=q 

We set uj(m) := m + and we denote by Ti(t), 1 < i < uj(m), the components of Y(t) 

labelled so that Yi(t) = Ui(t), VI < i < m. For u G and 1 < i,j < uj(m) we define the 
covariances 

Tijiu) := E[Yi{0)Y,{u)]. 

Observe that there exists a positive constant K such that 

I I < K'ild{u), Vi,j = 1 ,... ,u}{m), u G M"*, 

where if: is the function defined in ( 2 . 2 ). 


(3.10) 
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Using the Diagram Formula (see e.g.[19, Cor. 5.5] or [16, Thm. 7.33]) we deduce that for 
any 7 , 7 ' G Jm such that [ 7 ] = [ 7 '] = q there exists a universal homogeneous polynomial of 
degree q, in the variables Tij{u) such that 

E[H^{Y{0))Hy{Y{u))] = Tij {u )). 

Hence 


r ^ / I I 

V,,n = {2N)-^ ^ a(7)a(7') / {Ti,{u)) H (l - 


7 , 7 'el: 
I7l=l7'l=g 


du. 


(3.11) 


= :Rn{i,Y) 


From (3.10) we deduce that for any 7 , 7 ' G 3m such that [ 7 ] = [ 7 '] = q there exists a 
constant > 0 such that 

I P.^,y(Fij(rt)) I < C^^yilj{uy, Vit G M™. 

Since i/j G LP{MP), Vp > 1, we deduce from the dominated convergence theorem that 


lim i?^( 7 , 7 ') = i?^( 7 , 7 ') := f P^^y {rij{u))du, (3.12) 

N^OO 

and thus 

lim Vq^N = yq,oo -■= V' a{'j)a{j')R^{y,y') = [ E[pq{0)pq{u)]du. (3.13) 

N^OO ^ J^m 

7,7 SJm 

h\=h'\=<i 

Since Vq^N > 0, Vg, N, we have 

Vg,oo > 0 , Vq'. 


Lemma 3.4. For any positive integer Q we set 

y>Q,N ■■= Vq^N- 

q>Q 

Then 

lim (supU>Q,Ar) =0, (3.14) 

Q->oo jq 

the series 

^9,00 

9>1 

is convergent and, if 14 o is its sum, then 

Uoo = lim Vn = lim ^Vq^N- (3.15) 

N^oo N^oo 

9>1 

Proof. For s G we denote by Og the shift operator associated with the field X, i.e., 

0 ,X(.) =X(. + s). 

This extends to a unitary map L?‘{Tl) —)• L^(D) that commutes with the chaos decomposition 
of L^(D). Moreover, for any parallelipiped T we have 

Z{T + s) = 9sZ{T). 

If we denote by the set of lattice points 


:= [-iV, iV)”* n Z™ 
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then we deduce 

C(C^) = (2iV)-™/2 ^ esC{B), B = [0,1]"^. 

S^J^rn 

We denote by J’>q the projection 

^><9 ~ X!/ 

q>Q 

where we recall that IPg denotes the projection on the g-th chaos component of L‘^[Q). We 
have 

y>QCiCN) = {2N)-^/^ Y. (^sy>QC{B). 

S^Cm 

Using the stationarity of X we deduce 

V>Q^N = e[\T^qC{Cn)\^] =(2iV)— i.{s,N)E[?^q({B) ■ 0sT>qC{B)], (3.16) 

where z^(s, A^) denotes the number of lattice points t G Ljsf such that t — s G Ln. Clearly 

v{s,N) <{2N)^. (3.17) 

With K denoting the positive constant in (3.10) we choose positive numbers a,p such that 

'tp{s) < p< V|s|oo > a. 

K 

We split into two parts 

V>Q,N = y>Q,N + ^>Q,N^ 

where n i® niade up of the terms in (3.16) corresponding to lattice points s S L 2 N such 
that |s|oo < a + 1, while V”q jy corresponds to lattice points s £ L 2 N such that |s|oo > 0 + 1. 
We deduce from (3.17) that for 2M > a + 1 we have 

VYj^\ < i2N)-^{2a + 2r{2NrE[\?^QaB)f 

As Q —>■ 00 , the right-hand side of the above inequality goes to 0 uniformly with respect to 

N. 

To estimate V”q ^ observe that for s £ L 2 N such that |s|oo > a + 1 we have 

E[yyQC{B) ■ 9sy>QCiB)] = Y L iE[pg{t)]p,iu + s)]dtdu, (3.18) 

q>Q' 

where we recall from (3.9) that 


B JB 


p,{t)= Y := X v{m) 

7e3m, \'l\=q 


m(m + 1) 


Thus 


'y&m, \'y\=q 


E[pq{t)pqiu + S)] = E ( Y “(7)^^' 

763m, |7l=<? 

Arcones’ inequality [4, Lemma 1] implies that 

E[pq{t)pq{u + s)] <K'^'il;{s + u-ty Y I“(7)p7! 


763m, \l\=q 


(3.19) 
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The series divergent because the series Ylj&m ^ = {U,B), is 

the Hermite series decomposition of the distribution 6o{^/d^U)\ det A|. 

On the other hand, for 7 = {a, (3) G Jm we have 0 ( 7 ) = d{a)c{j3), where, according to 
(3.7) we have d{a) = ■^{2Trdm)~~Ha{0). Recalling that 


772.(0) = 772.+i(0) = 0. 


we deduce that 


and 


( 20 ! ^^ 2 .( 0 ) 


2 _ 1 / 2 r 


< 1 , 


d{afa\ <C = 


(27rdm)”*/^ 


This allows us to conclude that 


|a(7)p7! < {27Tdm)-^^\^ 0/3)'/3! < [ | det . 

763 ^, h\=q /36No^’"\|/3|<q 

Using this in (3.18) and (3.19) we deduce 

E[7^QC{B)-es‘P>QC{B)] 


Hence 


< CE„hr, 

t>Tr?. 


4|det^|2] [ [ ^{ 


s + u — tydudt 


Kq.n \< c'(X. <rK’'p‘‘-^ 

q>Q 


E 




B JB 


'ip{s + u — t)dud1^ , 


where we have used the fact that for |s|oo > a + 1 , |w|, |7| < 1 we have 'ijj{s + u — t) < p. 
Since p < ^, the sum 

q>Q 

is the tail of a convergent power series. On the other hand, 


sG-^2 


|s|oc 


> 0 .+! 


II Ipis + u- t)dudt < Y J 

^ ^ se£ 2 n 


[- 1 , 1 ]^ 


■■I 


'ip{s + u) <2 / 'il>{u)du < 00 . 


This proves that supjv|U>Qjvl 8 °®® as Q —>• 00 and completes the proof of (3.14). 

The claim (3.15) follows immediately from (3.14). □ 


Lemma 3.5. The asymptotic variance Uoo ’i-s positive. More precisely, 

^ 2,00 > 0 . 

Proof. From (3.13) we deduce 

V 2 ,oo = / E[p2{0)p2{u)]du, 


(3.20) 


(3.21) 
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where, according to (3.8) we have 

P2{t)= d{a)H^{U{t))fn{A{t)). 

aGNg^,nENo 

|a|+n=2 

The second chaos decomposition f 2 {A) is a linear combination of the polynomials p{A) and 
q{A) dehned in (3.4) where v = hm- 

In the above sum the only nontrivial terms correspond to a = 0 or a = {26ii, 26i2, ..., 25im)-, 
i = 1,... ,m. In each of these latter cases we deduce from (3.6) that 

d{a) := -^d(O), d(0) {2tt djn)~ ^, 

and we conclude that 

P2{t) = h(0) [f2{A{t)) - ^ f; H2{ Uiit) )) 

i=l 


= d{0) {xp{A{t)) + yq{A{t)) - ^2( Ui{t ))). 

For uniformity we set 


2=1 


2 : := —- 


k{A) 


so that 


P 2 {t) = d(0)(" xp{A) + yq{A) + H 2 {Ui) ). 


2 = 1 


We first express the polynomials p{A) and q{A) in terms of Hermite polynomials. We set 

( 1 


CLij ■ — 


/ 3 h, 

1 


-CLii, I — J, 




o-iji i 7^ J- 


We have 


(tr A)^ = 3hm [ I = ‘ihm YMi + QhmY/ 


i<j 


+ l) + 6hmY^i(ki)Hi{djj). 


i i<j 

p{A) = (tr^)^ - m{m + 2)hm = Ki I + 6 ^ ^fi(aM)-ffi(ajj) - m{m - 1) 

\ * i<j 

We have 

trA^ = ShmY'^li + ‘2hmY^ij = 3hm^i?2(an) + 2hmY dd2{dij) +fn{'^ + ‘^)k 

i i<j i i<j 


q{A) = hm ( sY dd 2 idii) + 2Y H^idij) 

i i<j 


Fo{t) = YH2{Uiit)), Flit) = YM^iiik), F2it) = YH2{dijit)) 

i i i<j 


Dehne 
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F3(t) = ^ii'i(aii(t))i/i(%(0). 


i<j 


Thus 


P2{t) = d{0)[xhm (3Fi(t) + 6 F 3 (t) - m{m - 1)) + yhm (3Fi(t) + 2F2{t )) + zFo{t) 

E[Fi{0)] = E[F2{0)] = 0, 


m(m — 1 ) , ,, m(m — 1 ) 

= -£^[-^ 3 ( 0 )] =- - --E[aii ( 0 ) 022 ( 0 )] = 


6 


We set 
Then 


F^it) = F^it) - EiFsit)]. 


E[Fo{0)Fs{t)] = E[FoiO)Fs{t)], E[Fi{0)F3{t)] = E[Fi{0)F3{t)], 
E[F2{0)F3{t)] = E[F2{0)F3{t)], 

P 2 {t) = d(0)(3xhm {Fi{t) + 2 F 3 (t))) +yhm (3-Fi(t) + 2F2{t)) +zFo{t) 


(3.22) 

-" '-^ 

=:Xl(t) =.Z2(t) 

To estimate E^p2{0)p2{u)~^ we will rely on the following consequences of the Diagram For¬ 
mula [ 8 ], [16, Thm. 3.12], 

Lemma 3.6. Suppose that Xi, X 2 , X 3 , X 4 are centered Gaussian random variables such that 
E[Xf] = 1, E[XiXj] = dj, \/i,j = 1, 2,4, i ^ j. 

Then 

E[H4{X4)H4{X^)] =C12, 

E[H2{Xi)H2{X2)] =2c?2, E[H2[X4]H4{X2)H4{X3)] = 2ci2Ci3, 
E[Hi{X 4 )H 4 {X 2 )H 4 {X 3 )H 4 {X 4 )] = C12C34 + C13C24 + C14C23. 


(3.23a) 

(3.23b) 

(3.23c) 

□ 


To compute the expectations involved in (3.21) we need to know the covariances between 
aij(O) and djk{t). These are determined by the covariance kernel 

C{t) = T[p{X)], 

where T denotes the Fourier transform. For any ii,..., G {1,... m} we set 

We have 


1 


Um = ^duX{t). «n(i)= 

y dm yS/i' 

Recalling that the spectral measure p{dX) has the form 

p{X) = u;(|A|)dA, 

we introduce the functions 

(A) := • • • Ajj,t(;(|A|) 

and denote by their Fourier transforms 


^ dlX(t), dij{t) = ^=dlt.X{t). 


'\/hm 




*1 


...,,(t) = (271)—/2 f e-(*’^)M,,..,,(A)dA. 
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Then 


E[umu^{t)] = -^Cij = 


E[Ui{0)ajj{t)]= -E[ajj{0)Ui{t)] = - 


1 


's/^hryidy 


dr)_ 
zCijjit) = 






E[Ui{0)djk{t)]= -E[ajkiO)Ui{t)] = - 


's/h-mdr 


-Cijkit) 


-3'ijk{t) j 7^ k, 


y/hm dr 

E[aii{0)ajj{t)] = — Ciijjit) = -- 3^iijj{t), 
Ofim '^•^m 

[^ii( 0 )flj/c(^)] — 

V ‘J 






We have 


= — 3 ~ji i 7^ J, k ^ i. 

f Lm 


E[Fo{0)Fo{t)] = ^E[H2{Ui{0)H2{Ujm = 


(3.24a) 

(3.24b) 

(3.24c) 

(3.24d) 

(3.24e) 

(3.24f) 


Using the fact that the Fourier transform is an isometry and the equality Mfj = 
we deduce 


MiiMjj 


[ Jijit)‘^dt= [ Mij{\fdX = [ Mii{X)Mjj{X)dX 

Jr™ Jr™ Jr™ 


and thus, 


E[Fo{0)Foit)]dt = { ^ 


dn 


J3 


(3.25) 


L 2 


F;[Fo(0)Ti(t)] = Y,E[H2mO))H2{d,jm 


3hm.d 




Since TfF = MaMjjjj we deduce 




F[Fo(0)Fi(t)]dt = 


3hmdn 


Eii^’ 








3/i™, 

1 / L 2 


(3.26) 


Arguing similary we deduce 


/ F[Fo(0)Fi(t)]dt= / F[Fo(t)Fi(0)]dt. 

Jr™ Jr™ 

F[Fo(0)F2(t)] = F2(C/,(0))F2(%fc(t)) = 2 ^ ^^g-,^.,(t), 

,• V UmU-m 


i,j<k 

V2 


i,j<k 


[ F[Fo(0)F2(t)]dt=(^j;M,,,^j;M 
Jr™ \ , 


jjfcfc y 

j<k / 


F[Fo(t)F2(0)]dt. (3.27) 
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1 


£;[Fo(0)F3(t)] = E[H2{UimHi{ajj{t)Hi{akk{t))] = -2 ^ 

i,j<k i,J<k 

£:|F„(0)F3(()I* = E {Mijj, Mikk)]^2 


\n 


V 2 




^jjkk 


j<k 


3hmd n 


L 2 


i,j<k 

f E[Fo{t)F3{0)]dt. 


(3.28) 


We have 


E[F^{0)F,{t)] =Y,E[H2{du{0))H2{d,,{t))] = 






E [ Fi (0)Fi {t)]dt = {^Y ^ E 


3hr. 


(3.29) 


E[Fi{0)F2{t)] = Y E[H2{dumH2{djkit))] =^Y 

™- i,j<k 


i,j<k 


[ E[Fi{0)F2{t)]dt = ( Y^ Man, E Mjkjk 
JR^ \onm . j<k / 


V2 


(3.30) 


L 2 


i;[Fi(0)F3(t)] = Y E[H2{aii{Q))Hi{ajj{t))Hi{akk{i)] = {t)Jiikk{t). 

Ky / 

i,j<k ^ iij<k 


[ E[F, {0)Fs {t)]dt=(^Y ^ E ^Fkk ) • 

Jk- yilm . j^k 7^2 

£;[F2(0)F2(t)] = E[H2{d.,mH2{dk,{t)]=2 Y 


(3.31) 


i<j,k<£ 


i<j,k<£ 


(3.32) 


/ £;[F2(0)F2(t)]dt = (pj;M,,,,-,pj;M,,,,) . 

\ Kj i<J / 12 

F[F3(0)f3(t)] = E F[if2(Si3(0))ffife(())ff3(a„(())] =-|- 3~ ijkk{^)3'ij£i{t)^ 

i<j,k<i ™ i<j,k<i 


[ E[F2{0)Fs{t)]=(^Y^^ 

JR"» \ , . 


V 2 


oL Yl^kkil ) 
fc<£ / 


(3.33) 


£;[F3(0)F3(t)] = ^ £;[2f?i(aii(0))2H'i(%(0))2F/i(afcfc(t))2f?i(S,,(t))] 

i<j^k<i 

= ^ £^[aM(0)ajj(0)]£^[afcfc(t)a^£(t)] + ^ (3'iifcfc(i)3j,7T^(i) + 3'M££(i)3jjfcfc(i)) 

i<j^k<i ^ i<j^k<i 

= E[Fs{0)f + ^ Y ( 3 "iikk (^) 3 ~ jjiiij') + 3^u£e{t)3^jjkkit)^ ■ 

i<j,k<l 

E[FmF^{t)] = E[FmFz{t)]-E[Fm? = ^ E • 


i<j,k<£ 
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We deduce 


f E[F3{0)F3{t)]dt =^ [ (3'iikkit)3'jjuit)+3'nuit)3'jjkkit))dt 




and we conclude 


E[F3{0)F3{t)]dt =^ E 


(3.34) 


i<j ™ k<e / ^2 


To put the above equalities in perspective we introduce the functions 

= = G2 = ^Y^^Fkk, G3 = \g2. 


7 / —to J. OT / —7 / ^ o' 

i i j<k 

The assumption A3 implies that Go,Gi,G 2 S Using the notation 

Fi.Fj:= [ E[FM,Fjit)]. 

JR^ 

we can rewrite the equalities (3.25, ..., 3.34) in a more concise form 

Fi • Fj = Fj • Fi = {Gi, Gj)L 2 , Fi» Fs = Fi» Fs, Vi, j = 0,1,2, 

F 3 • Fi = Fi • F 3 = {Gi, G 3 )l 2 , Vi = 0,..., 3. 

From (3.21) and (3.22) we deduce 

V2,oo = f P2{t)dt = d{0)'^ (sxhmZi + yhmZ2 + zFo)'\ • (sxhmZi + yhmZ2 + zFo)'] 

liam V / V / 


— d(0)^ 3xhjyi{Gi + 2 G 3 ) + yh^i^Gi + 2 G 2 ) + 2;Go 


L 2 


((i( 0 ))^ 3xhm + 3yhM ^Gi + ^^ 2 ^ + 


2 

L 2 


— (i(0)^ 3hm{x + y) + 2 ^* 2 ^ + zGq 


2 

L 2 ' 


The functions Gi + |G 2 and Go are linearly independent and 

^ = -^MA) = -iT;[|detA|] / 0 . 

Hence V 2 ,cxd > 0. 

This concludes the proof of Proposition 3.3. 


□ 

□ 


Remark 3.7. The numbers x,y that describe f 2 {A), the 2nd chaos component of f{A) seem 
hard to compute in general. In Appendix A we describe their large m asymptotics; see (A. 16). 

□ 
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3.4. Conclusion. As explained in Subsection 2.3, to conclude the proof of Theorem 2.2 it 
suffices to establish the asymptotic normality of the sequences 

This follows from the fourth-moment theorem [26, Thm. 5.2.7], [29]. Here are the details. 

Recall from [16, IV. 1] that we have a surjective isometry 0^ : , where 

is the g-th symmetric power and is the g-th chaos component of . The multiple Ito 
integral Iq is then the map 

We can write QN,q as a multiple Ito integral 

Ca',, = L[<], 

According to [26, Thm.5.2.7(v)], to prove that CN,q converge in law to a normal variable it 
suffices to show that 


lim llfff ( 8 )r S-f ll.r®( 2 i- 2 -) = 0> Vr = 1,... - 1. 


N^oo 


(3.35) 


To show this we invoke arguments following the inequality (18) in the proof of [13, Prop. 2.4] 
which extend without any modifications to the setup in our paper. 

Appendix A. Asymptotics of some Gaussian integrals 

We want to give an approximate description of the 2nd chaos component of f{A) = | det A\ 
when m ^ 0 . 

Observe that if u : —)• K is a continuous function, homogeneous of degree k, then for 

any u > 0 we have 

Ei,-Ju{A)] = {2v)^E i/2[u{A)]. 


Proposition A.l. Set Cm ■= 22r We have the following asymptotic estimates as 

2 


m 


oc 


E^y2[f{A)] r^Cm^m i 


E^y.[p{A)f{A)] 
Es[q{A)f{A)] r 


TT 

2 Cm 


3 

m 2 . 


Cr 




7 

m 2 . 


(A.l) 

(A.2) 

(A.3) 


Proof. We need to make a brief detour in the world of random matrices. 

We have a Weyl integration formula [3] which states that if / : Sm —>• K is a measurable 
function which is invariant under conjugation, then the value f{A) at A G §>m depends only 
on the eigenvalues Ai(A) < • • • < An(A) of A and we have 


-E^goe^(/(V)) = 


1 


Zm{v) 


/(Ai,...,A^) ( H |V-A,| I ne-^|hAi---dA^|, 


2=1 


— ■Qm,v (•^) 


(A.4) 
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where Zm{v) can be computed via Selberg integrals, [3, Eq. (2.5.11)], and we have 


= (2„)=^Z„, Z„ = (2^)fm!n^=2?m!nr(:^). (A,5) 

j.lU2) j.l V2/ 

For any positive integer n we define the normalized 1-point correlation function pn,v{x) of 
GOE;; to be 

Pn,v{x') — 7 r / Qn,v{x^ A 2 , • • • , ^n)dXi • • • dXn- 

Zn[V) 

For any Borel measurable function / : M —)• M we have [12, §4.4] 

-E;GOE”(tr/(X)) = / fiX)pnA^)dX. (A.6) 

^ Jr 

The equality (A. 6 ) characterizes pn,v Let us observe that for any constant c > 0, if 

A e GOE); ^cA e GOE^. 

Hence for any Borel set H C M we have 

/ Pn,Av{x)dx= / pn,v{y)dy, 

JcB Jb 

and we conclude that 

cPn,c^vicy) = Pn,v{y)^ Vn, c, y. (A.7) 

The behavior of pn,v for large n is described the the celebrated Wigner semicircle theorem. 

Theorem A.2 (Wigner). As n 00 , the probability measures 

Pn,vn-^iX)\dX\ = n^^^pn,v{n^^^X)\dX\ 

converge weakly to the semicircle distribution 

poo,vW\dX\ = I\x\<2^^V^v-^\dX\. □ 

We have the following result of Y. Fyodorov [14]; see also [22, Lemmas C.l, C2.]. 

Lemma A.3. Suppose v > 0. Then for any A G M we have 

m+l 


r(i) 


-E'GOe;;, (I det(A-I-Hj) = (2u) 2 Cme-i^Pm+i,v{X), Cm-=‘^^T 


E,.^{\det{A)\) = {2v)^ ^ 


V 2 


'Ey 


-E/GOEJi, (I det(A-|-H] )e ^vdX 


/„ . EL+i Cm 

= { 2 v} 2 


V2 


TTV 


J 

JR 


Pm+i,viX)e 4„dA. 


(A. 8 ) 


(A.9) 


□ 


We will also need the following asymptotic estimates. 
Lemma A.4. Let k be a nonnegative integer. Then 
1 




f 


2 fc „ v^V^{2k-l)\\ 

e 2 i 2 dX ~- 1 = - as n ^ 00 


(A.IO) 














A CENTRAL LIMIT THEOREM 


21 


Proof. Consider the function 


Then 


<,2k 


w{X) = 


e 




I 

JR 


w{X)dX = 1. 


We set Wn{X) := y/nw(y/nX). The probability measures Wn{X)dX converge to Jo and we have 


2k + l 


p^i{X)X^'^e-^dX = 


V 2 (2k — 1 )!! 


Jr ”’2 Vn 

Arguing exactly as in [23, Sec. 4.6] we deduce 


Pn,J-W'Wn{X)dX. 


lim / PnJ_{X)wn{X)dX = p^i{0) = J-. 
n^oo Jjg ’2n °°’2 V TT 


□ 


The estimate (A.l) follows from (A.9) and (A.10). 

To simplify the notation we set 

Es-.= E^ 1 J 2 . 

Let us observe that the equality (2.5) implies that 


^ / F;G[p(A + i?)/(A + i?)]e-^'dA, 

Jr 

(A.lla) 

^ / F;G[g(A + i?)/(A + i?)]e-^'dA. 

JR 

(A.llb) 


Es[q{A)f{A)] = 

To estimate Es\p{A)f{A) ] and Es [ q{A)f{A) ] for m large we use a nice trick we learned 
from Yan Fyodorov. Introduce the functions 

4 ), T : M X (- 00 , 1 ) M, 

$(A,z) := T;G[|det(A + S)|e^*’’(^+^)], T(A, 2 ) ;= FIg [ | det(A + 5)|ei ]. 

Obviously 

$(A, 0 ) = T(A, 0) = FIg [ I det(A + 5)| ] 
so both T(A,0) and T(A,0) can be determined using (A. 8 ). 

Observe next that 

ch",(A,0) = F;G[(tr(A + i?))'|det(A + i?)|] = Eg[p{X +B)f{X +B)], (A.12a) 

2T",(A,0) =F;G[(tr(A + B)2|det(A + S)|] = FIg [ ^(A + B)/(A + S) ]. (A.12b) 

We have the following key observation. 

Lemma A.5 (Y. Fyodorov). 

<I)(A, z) = + z, 0) = e ~^+ z, 0) 


T(A,^) = 


e 2 (i-^) 


m(m+3) 

{I- Z) 4 




VT^ 


:,0 


(A. 13a) 
(A.13b) 
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Proof. Using (2.4) we deduce 

^X,z) = Km [ + = f e"^ 

J Sm 

(make the change in variables C := B — z) 

= [ I det(A + z + dB 

J Sm 

= [ I det(A + z + C)\] = + ^, o). 

Similarly, for z < 1 we have 

'^{X,z) = Km [ \det{X +B)\e^ 2 ^^(^+B)^-i^^B^dB 


mz\^ f , 

= 62 / I det(A + B)\e 2 

Making the change in variables B = (1 — z)~^^'^C so that 


'dB. 


dB = {1 — z) 4 dC det(A + i?) = (1 — u) 2 det(A\/l — 2 + C). 


We deduce 


(c-- 


\z 


4/(A,^) = 


B) 


6 2 


Kr 


7Tl(7Tl + 3) ^ 

(1 — Z) 4 


det(A\/l — z + C)\e ^ 




dC 


e 2 


m(m+3) 

(1 — Z) 4 


det(A\/l — z -\— + B)\e 2 ^^^^dB 

\/\ — z 


6 2 


m(m+3) 

( 1 - 2 :) 4 


qj 


Vi — 


□ 


The asymptotic behavior of Es[p{A)f{A)~\. Using (A.13a) we deduce 

^.L(A,0) = e-=^^C|„„(ef(^«)=^.(A + .,0))=e-=^^|,(e=^^.(A,0)). 


Using (A. 11a) and (A. 12a) we deduce 


Bs[p{A)f{A)]=-^J^e "*2 ^(e ’"2 $(A,0))e ^^dX 

1 f Cp' / mA^ m + 2 \2 

2 ( , ))e 2 

Since <b(A, 0) has polynomial growth in A, we can integrate by parts in the above equality 
and we deduce 

Es[p{A)f{A)] = ^^e^$(A,0)^(e"^^')o!A 
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(A.8) Cm f 


f mA^ ^ /,x U / rn+2 ^2 x „ 

■ )<iA 


Cm r (yn,+ l)A d . _m+2;^2 A 

-V^L° ‘ '’■"+'3Wsp(' “ 

/ Pm+i4(A)((™ + 2)A^-l)e-AdA (0), 


We have thus proved that 


y/m + 1 


Es[p{A)f{A)] ~ + 2)("i + l)/5oo,i(0) asm-^oo. (A.14) 


This proves (A.2). 

The asymptotic behavior of £^5 [ g(A)/(A) ]. We set 


m(z) : = 


e2(i-2) 


m(m+3) 

(1 -^) 4 


mA 2 

= e 2 - 


m(m+3) 

(1 - Z) 4 


Then 


T(A, z) = u(z)T(A(l - 0), 

t;(A,z) = u'(z)T(A(l - + ^u{z){l - z)-^/H'^{X{l - z)-^/\Q). 

Kz{\z) = u"(^)T(A(l - z)-^/^0) + ^v!{z){l - z)-3/2^r'^(A(l - z)-^/^,Q) 

+ ^^(«(^)(1 - )'hUA(l - + ^-^u{z){l - z)-%';,(A(l - ^)-'/^0). 


Thus 


A 


'h;',(A,0) = u"(0)'h(A,0) + -n'(0)T',(A,0) 
A 8 8A^ 

A 


— u"(0)'I'(A, 0 ) + ^ ( 2m'(0) + 2^(0) ) '^a(A, 0) H (A, O). 

Setting K{m) = we deduce 

u'{z) = f e^(l - 


dz \ 


mA 2 / mX^ 
= e 2 


^e2(Tvi) (1 _ ^ ^(j„)e2(TVi) (1 _ j . 


Thus 


We set 


,, , mX^ , , 
n (0) = — -h K(m). 

X / mX^ 


-Ai (A) = - (n (0) + -) = - ( ^ + «:(m) + - ). 


Similarly, we deduce 


«"(0 ) = ^ + k("i) + 2^ + K(m) ^ + K(m) + 1) 
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m 


2^4 


+ {K{m) + l)m\^ + K{m) ( K{m) + 1 ) 


= :44lo(A) 


We set A 2 (A) := We have 

0) = A 2 (A)^';(;,(A, 0) + Ali(A)'I/'^(A, 0) + Alo(A)'I/(A, 0). 

Using (A.lib) and (A.12b) we deduce 

Es [ q{A)f{A) ] = ^ ^ ( Al 2 (A)'h'^A(A, 0) + Ai(A)^';,(A, 0) + Ao(A)^(A, 0)) e'^'dA 


1 




= g= / 4,(A,0))|^(A2(A)e-"’)-A(A,(A)e-"’) + A„(A)e-"“ ) dA. 


vr 


dA 2 
1 


P4,m(A)$(A)^'h(A,0)e-^ dA, 


where 


-P4,m(A) — A2(A) — 4AA2(A) + 4A^A2(A) — A]^(A) + 2AAi(A) + Ao(A) 

= C' 4 (m)A^ + C 2 {m)X^ + Co{m), 

where the coefficients CQ{'m),C 2 {m),C 4 {m) are polynomials in m. Recalling that 

'I'(A,0) = Eg[ I det(A + B)\ ]C'mPm+i,i/ 2 (A), 

we deduce 


Es[q{A)f{A)] = 


Ciirn) f M^\4e , C2{m) I 




Pm+i,i(A)A /^Pm+i,i(A)A 


TT y/-K 

C'o(m) f 


Using (A. 10) with u = 1/2 we deduce that as m —)■ oo we have 




Upon investigating the definition of Ao(A), Ai(A), and Ao(A) we see that of the three 

deg Co(m) = 4 > deg C 2 (m), deg C 4 (m). 

The degree-4 term in Co(m) comes from the product 

4 

Tin 

2 K(m)(«;(m) -|- 1) = —h lower order terms. 

We conclude that as m —)• oo we have 

Es[q{A)f{A)] ~ mi 

V ZTT 

□ 
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To understand the 2nd chaos component of | det A\ we need to also understand the inner 
product in For simplicity will write E instead of the more precise Eg-i^. 

We know that 

E[p{A)]= [g(A)] = m{m + 2)v. 

This implies that 

E[p{A)'^~\ = E[{tr A)^~\= + 2 )‘^v‘^. 

To compute E[p{A)q{A)], E[q{A)'^] we will use Wick’s formula, [16, Thm. 1.28]. We have 


p{A)q{A) = ['^al + 2'^, 


i<j 


\ k k<i 


Oh 


X] ^ X/^ X/ + X/ + 4 ^ auajjah. 


i<k 


i, k<i 


k, i<j 


i<j, k<£ 


Si 


S 2 


We have 


E[S3] = E 
T^[55] = E 


T;[5i] = E 


2 Y1 

i, k<i 


S3 




54 


5s 


= mE[a\]\ = 27mv‘^. 


4 anajjah 

i<j, k<i 


= w?{m — l)E[a\i]E[a\2] = 3m‘^{m — l)i;^. 

= nn?{m — l)^Tl[aiia 22 ]T^[ai 2 ] = nn?{nn — l)^u^. 


E[S2] = E 


2 


i<k 


= m{m - l)T;[afia^ 2 ] 


Using Wick’s formula we deduce 

-^[^ 11 ^ 22 ] ~ -^[®ii]-^[®22] + 2£l[aiia22]^ = llu^. 

Hence 

T1[52] = llm(m — l)u^. 


(A.15) 


E[Si] = E 


2 'y ^ aj^i^auajj 

k, i<j 


2E 


+ 2E 


+ 2E 

yy ^kk^ii^jj 


_i<j 


J<j 


i<j, k^i,j 


= 4E 


E 

Kj 




+ m{m — l)(m — 2 )Tl[aiia 22 a 33 ] 


= 2m{m — l)Fl[a^^a 22 ] + +m{m — l)(m — 2 )Tl[aiia 22033 ]- 
Using Wick’s formula we deduce 

Tl[afia22] = 3£'[a3]^]Tl[aiia22] = 

£^[011022033] = £[aii 022 ]£[a 33 ] + 2£[oiia22]^ = 5 u^, 
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Hence 


E[S 4 ] = 18m{m — l)v'^ + 5m(m — l)(m — 2)v'^ = m{m — l)(5m + 8)v^. 


We have 


/ 


q{A)^ = 


\ 


i k<i 


\ 


= X^ + Y^ + 2XY. 




X Y 

The random variables X and Y are independent and thus 

E[q{A)'^] =E[X'^] + ElY"^] + 2E[X]E[Y]. 

We have 

E[X] = 3mv, E[Y] = m{m — l)^, 2E[XY] = 6m^(m — l)v^. 

Next, 

^ afi + 2 ^ aW 




Kj 


(A IS") 

= mE[a\^ + m{m — l)E[a‘fia 22 ] = 27mv‘^ + llm(m — + 16mu^, 


r^ = 4W 






2 2 
j ^kl 


Kk<i 


k<l 


i<j, k<i 


m 


E[Y^] = A[ ^ U;[at2] + 8 


(T) 

2 


E[a 


2 12 
I 2 J • 


= 6m(m — ^ ~ m{m — l)u^ (6 + 2(m + l)(m — 2)). 

We summarize the results we have obtained so far. Below we denote by o(l) a function of m, 
independent of v that goes to 0 as m —)• 00 . 

£;[p(H)] = m{m + 2)u, 

] = 3w?{m + 2)‘^v^ = 3m^u^(l + o(l)), 

E[q{A)] = m{m + 2)u, 

= mv^{2m^ + 2w? + 9m + 14) = 2m'^u^(l + o(l)), 

E\p{A)q{A)] = (m^ + 3m^ + 12m + = m^u^(l + o(l)). 

We have 

E\p{A)‘^] = E\p{A)‘^] — 2m(m + 2)Tl[p(H)] + m(m + 2)u = 2m^u^(l + o(l)). 
E[p{A)q{A) ] = E\p{A)q{A)] — m^(m + 2)^u^ = —m?v^{l + o(l)), 

E[q{A)‘^\ = mv‘^v^{l + o(l)). 

Thus, in the basis p{A), q{A) of the inner product is given by the symmetric matrix 

2 0(1) ^ 


4„,2 


‘2” = ”” [„(!) 1 

This proves that the component of f{A) in has a decomposition 

/2(2l) = Xmp{A) + ymqmq{A), 
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where, as m —>■ oo 


1 


Xm. 


2m^v'^ 


/ \ '^+2 
{2v) 2 

2m^v‘^ 

1 


Es-^ [p{A)f{A) ] - m{m + 2)vEgv^ [ f{A) ] 

E^,nlp(A)f(A)]-^A^E,.JnA)]). 


ym 


m^v‘^ 


Es-^[piA)f{A)] -mim + 2)vEsv^[f{A)] ), 


/ \ m+2 

{2v) 2 
2m‘^v'^ 


E ^/ 2 [q{A)f{A)] - 


m{m + 2) 


Bs:.[/0)] 


Using (A.1),(A.2) and (A.3) we deduce that there exist two universal constant zi, Z 2 ,independent 
of m and v such that, as m —)■ oo, 


ziCmV 2 rn 


^rr,-5/2 


Vm 


Z 2 CmV 2 m 


(A.16) 


In the problem investigated in this paper the variance v also depends on m, u = hm- Recall 
that the constant Cm grows really fast as m —)• oo 

log Cm ~ ^mlogm. 

Proposition A.6. The Gaussian vector 

J2{X) := (A(0),VA(0),V2a(0)). 

is nondegenerate. 


Proof. We set H := V^(0) and we denote by Hij its entries. The equality (2.8d) shows that 
H G Is a centered Gaussian random real symmetric matrix whose statistic is dehned by 
the equalities 

E[E^,]=3hm, E ]=E[Hfj]= hm, Vi / j, 

while all the other covariances are trivial. This shows that the second jet J 2 (X) is the direct 
sum of mutually independent Gaussian vectors, J 2 {X) = A 0 Hq 0 D, where D = VA(0), 
Hq is the vector with independent entries and A is the vector 

A={X{0),Hu,...,Hmm). 

The components Hq and D are obviously nondegenerate Gaussian vectors. Thus, the jet 
J 2 {X) is nondegenerate if and only if the component A is. The covariance matrix of A is 
Rm{sm, dm, hm) where for any s,d,h > 0 we denote by Rm{s, d, h) symmetric (m+l) x {m+1) 
matrix with entries 


roo = s, roi = -d, Vi = 1,.. ., m, ru = 3h, rij = h, VI < i < j < m. 

Note that multiplying the hrst row by and then the hrst column by we deduce 

det Rm{s, h,d) = s det Rm{i,d, h), d = ds~^^'^. 

If we add the first column multiplied by d to the other columns we deduce that 

det d, h) = det Gm{^h — (P,h — (P), 

where Gm{(^,b) denotes the symmetric mx m matrix whose diagonal entries are equal to a, 
and the off diagonal entries equal to b. As explained in [23, Appendix B], we have 

det Gm{a) = (o — (a 0 (m — 1)6). 
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Thus 

det Rmis, d, h) = s{2h)'^~^ ((m + 2)h — mcP ) = ((m + 2)hs — md? ). 

Thus J 2 {X) is nondegenerate if and only if Using (2.7) we deduce that 

dmSm _ ^ 

m + 2 Im+iiw)^ 

From the Cauchy inequality we deduce that Im+i{w)‘^ < Im-i{w)Im+ 3 {w)- We cannot have 
equality because the functions ^w{r) r^~ and w{r) r“ 2 “ are linearly independent. □ 
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